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SECTION  1 
INTRODUCTION 

Structural  sandwich  composites  form  a  very  important  part 
of  the  total  materials  pool  available  to  structural  designers. 

The  high  strength-to-weight  ratio  which  is  characteristic  of 
structural  composites  often  gives  them  a  favorable  position  over 
more  conventional  materials  for  structural  applications  for  which 
weight  is  a  limiting  or  controlling  factor  in  the  material 
selection.  Consequently,  structural  sandwich  composites  are  used 
widely  for  various  applications;  in  particular,  they  are 
especially  valuable  for  aerospace  applications  due  to  their  high 
strength  and  low  weight  traits. 

Since  structural  sandwich  composites  are  used  in  many 
structures  for  which  low  weight  is  a  desirable  property,  they 
are  natural  candidates  for  application  of  mathematical  optimum 
design  techniques.  The  consideration  of  sandwich  in  minimum 
weight  design  studies  offers  greater  challenge  and  wider 
opportunity  for  significant  weight  reduction,  because  additional 
design  parameters  are  available  for  adjustment  by  the  structural 
designer.  For  example,  in  addition  to  the  normal  design 
parameters  such  as  panel  thickness  and  material,  the  designer 
of  sandwich  panels  can  choose  the  thickness  of  individual  faces, 
the  thickness  of  the  core,  the  type  of  core,  the  core  material, 
as  well  as  parameters  which  control  the  density  of  the  core. 

Although  the  availability  of  numerous  design  variables 
offers  the  creative  designer  opportunities  which  would  not  normally 
exist,  the  process  of  optimum  design  of  sandwich  is  necessarily 
more  complex.  The  purpose  of  this  report  is  to  present  an 
approach  to  the  minimum  weight  design  of  sandwich  panels. 

Section  1  defines  the  scope  of  the  study,  and  presents  a  body 
of  literature  concerned  with  the  design  of  structural  sandwich 
composites.  Section  2  is  concerned  with  a  development  of  the 
mathematical  representation  of  the  problem  of  the  minimum 


weight  design  of  sandwich  panels.  The  general  constrained 
minimization  technique  used  is  discussed.  Also,  the  mathematical 
representation  of  the  sandwich  design  parameters,  weight,  and 
design  constraints  within  the  framework  of  the  selected 
optimization  procedure  are  presented.  A  description  of  the 
computer  program  which  was  developed  to  implement  the  sandwich 
panel  minimum  weight  design  procedure  is  presented  in  Section  3. 
Instructions  for  operation  of  the  interactive  program  are  given, 
and  typical  applications  of  the  program  to  representative  design 
examples  are  considered.  Section  4  concludes  that  the  results 
obtained  in  this  panel  optimization  study  are  very  encouraging, 
and  recommends  that  efforts  be  undertaken  to  apply  minimization 
techniques  to  more  general  structures  which  contain  sandwich. 

1 . 1  SCOPE 

This  report  is  concerned  with  the  automated  minimum 
weight  design  of  flat,  rectangular,  sandwich  panels  subjected 
to  edgewise  compression  loads  {Paragraph  2.3  gives  a  more 
complete  description  of  the  particular  geometrical,  material, 
and  loading  configurations  considered) .  The  rather  narrow 
problem  scope  was  selected  as  a  basis  for  demonstrating  that 
application  of  mathematical  constrained  optimization  techniques 
to  the  design  of  sandwich  construction  is  both  feasible  and 
practical.  The  sandwich  panel  design  problem  considered  contains 
all  of  the  basic  ingredients  of  a  more  complex  problem  for 
which  sandwich  panels  are  only  individual  components.  That  is, 
design  parameters  and  design  parameter  constraints  are  used 
which  are  not  normally  considered  in  more  conventional  designs. 

It  is  recommended  in  Section  4  that  further  work  is  warranted 
in  combining  the  advantages  of  mathematical  programming  approaches 
in  individual  panel  design  (this  study) ,  with  the  optimality  design 
techniques  used  for  overall  design  of  more  extensive  structures 
which  contain  panels  as  individual  components. 
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1.2  LITERATURE 


Quite  a  large  body  of  literature  exists  concerning  the 
general  topic  of  optimum  structural  design  and  the  more  specific 
subject  of  optimum  sandwich  composites  design.  A  complete 
review  of  the  pertinent  literature  is  not  attempted  here. 

Instead  we  merely  provide  a  bibliography  of  sources  (Appendix  A) 
which  pertain  to  the  subject.  The  bibliographic  list  has  been 
compiled  from  the  Structural  Sandwich  Composites  Bibliography 
of  Reference  1. 
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SECTION  2 

TECHNICAL  APPROACH 

This  section  presents  the  mathematical  basis  for  the 
computer  program  described  in  Section  3.  The  sandwich 
optimization  study  considered  in  this  report  is  put  into  the 
general  framework  of  a  general  constrained  optimization  problem 
(Paragraph  2.1),  and  solutions  are  obtained  using  an  approach 
called  the  sequential  unconstrained  minimization  technique 
(Paragraph  2.2).  The  particular  geometry,  material,  and  loading 
of  the  sandwich  panel  studied  are  defined  in  Paragraph  2.3,  and 
the  parameters  chosen  for  performing  optimization  studies  are 
defined  in  Paragraph  2.4.  The  equations  which  quantify  the 
objective  function  (weight  in  this  study)  and  the  constraints 
on  the  design  parameters  are  given  in  Paragraphs  2.5  and  2.6. 

2.1  GENERAL  CONSTRAINED  OPTIMIZATION  PROBLEM 

The  constrained  optimization  approach  has  been  selected 
for  this  sandwich  panel  optimization  study.  A  complete  account 
of  this  method  is  found  in  References  2  and  3.  a  brief  description 
of  constrained  optimization  is  presented  here  for  completeness. 

A  statement  of  a  constrained  optimization  problem  in 
mathematical  terms  is: 

Select  design  parameters  d]_,...,dn  such  that  an 
objective  function  W(d^,...,dn)  is  minimized,  subject 
to  the  inequality  design  parameter  constraint  functions,  (2.1) 
9i (di » • • • ,dn)  ^  0,  i  =  1 , . . . , N ,  and  the  equality  design 
parameter  constraints,  (d]_ , . . .  ,dn)  =  0,  i  =  1,...,M. 

In  physical  terms,  the  n  design  parameters  d^,...,dn  can 
be  viewed  as  the  axes  of  an  n-dimensional ,  orthogonal  coordinate 
space  on  which  the  objective  function  W(d^,...,dn)  is  defined. 

Each  point  of  this  n-dimensional  space  represents  a  possible  set 
of  design  parameters.  The  constraint  functions  serve  to  divide 
the  total  design  parameter  space  into  acceptable  and  unacceptable 


regions  as  far  as  valid  designs  are  concerned.  Then  the 
optimization  problem  can  be  visualized  as  the  selection  of  that 
point  in  the  acceptable  region  of  the  design  space  which 
corresponds  to  the  minimum  value  of  the  objective  function. 

There  are  essentially  two  general  approaches  which  can 
be  used  to  solve  the  optimization  problem  stated  above;  these 
basic  techniques  are  known  as  direct  methods  and  indirect 
methods.  The  direct  methods  employ  a  function  minimization 
scheme  which  operates  directly  on  the  objective  function,  with 
the  constraints  being  considered  as  limiting  surfaces.  Specialized 
techniques  have  been  developed  for  determining  whether  an 
optimum  has  been  attained  upon  encountering  constraints  and  for 
proceeding  with  the  optimization  if  necessary.  The  indirect 
methods,  on  the  other  hand,  rely  on  a  reformulation  of  the 
problem  which  converts  the  constrained  optimization  problem  into 
an  unconstrained  optimization  problem.  This  is  an  extremely 
attractive  approach  since  it  means  that  standard,  well  established 
unconstrained  minimization  techniques  can  be  used  to  obtain  an 
optimum  design. 

The  indirect  approach  has  been  chosen  for  the  optimization 
of  sandwich  panels  reported  here.  In  particular,  the  Sequential 
Unconstrained  Minimization  Technique  developed  by  Fiacco  and 
McCormick  has  been  used  to  obtain  the  set  of  design  parameters 
which  yield  an  optimum  weight  sandwich  component.  This  method 
is  described  briefly  in  the  following  paragraph. 

2.2  SEQUENTIAL  UNCONSTRAINED  MINIMIZATION  TECHNIQUE 

The  Sequential  Unconstrained  Minimization  Technique  (SUMT) 
has  been  chosen  due  to  its  proven  applicability  combined  with 
the  fact  that  it  is  simple  and  relatively  foolproof  to  use. 

The  mathematical  theory  behind  this  method  is  contained  in 
Reference  4  while  examples  of  the  application  of  the  technique 
to  various  optimization  problems  are  given  in  Reference  5.  The 
particular  form  of  the  SUMT  which  is  used  generates  intermediate 
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designs  which  all  lie  inside  the  acceptable  region  of  the  design 
space . 

The  basic  idea  behind  the  SUMT  is  relatively  simple.  In 
this  approach  the  objective  function  is  augmented  with  a  "penalty 
function"  which  is  designed  to  contain  the  effects  of  the  design 
constraints.  Then  a  sequence  of  unconstrained  minimizations  is 
performed  on  the  new  functions,  with  each  successive  minimization 
producing  a  result  which  is  closer  to  the  true  optimum. 

In  equation  form,  the  SUMT  objective  function  is 


$ (r ; d^ , . . . , dn )  =  W(d^, . . . ,dn) 


N 

r  E  In  g. (d. 
i=l  1  L 


'V 


+ 


M 

E 

i=l 


Hi  (di'*“  'V 


(2.2) 


where  r  is  a  parameter  which  controls  the  magnitude  of  the 
penalty  function  (the  last  two  terms  in  Equation  2.2);  if  r  =  0 
then  $  =  W.  The  general  idea  behind  the  SUMT  is  to  select  a 
value  for  r,  perform  a  minimization,  reduce  r,  perform  another 
minimization,  etc. ,  until  r  is  made  sufficiently  small  that 
min  $  *  min  W.  References  2 ,  4 ,  and  5  contain  more  complete 
discussions  of  the  theory  of  this  method. 

The  following  algorithm  has  been  used  in  applying  the 
Sequential  Unconstrained  Minimization  Technique  to  an 
optimization  problem  (refer  to  Figure  1) : 

(a)  Pick  a  starting  value  for  r  and  select  an  initial 
set  of  design  parameters  (dlf...,dn)0  satisfying 
all  the  constraints  (Equation  2.3). 

(b)  Minimize  <J>  (Equation  2.2)  to  obtain  (d. ,...,d  )„, 

i  n  m 

where  M  denotes  the  Mth  pass  through  the 
algorithm. 

(c)  Check  for  convergence  to  the  optimum. 
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Figure  1.  SUMT  Logic  Diagram. 


(d)  If  the  convergence  criterion  is  not  met,  reduce 
r  by  r  -4 —  cr ,  where  c  <  1 . 

(e)  Compute  a  new  starting  point  for  the  minimization, 
initialize  the  minimization  algorithm,  and  repeat 
from  Step  1. 


2.3  STRUCTURE  AND  LOADING 

The  class  of  structure  and  loading  considered  for  this 
initial  sandwich  optimization  study  is  that  represented  in 
Chapter  5  of  the  MIL-HDBK-23A. ®  In  particular,  the 
reported  study  is  directed  at  simply-supported,  rectangular, 
three  layer,  sandwich  plates  subjected  to  inplane  compression 
loads. 


The  geometry  of  the  subject  sandwich  plates  is  shown  in 
Figure  2.  The  planform  dimensions  are  denoted  by  a  and  b.  The 
two  face  plates  are  of  uniform  but  different  thicknesses  t^ 
and  t2 ,  and  the  core  has  uniform  thickness  tc.  The  total 
thickness  of  the  sandwich  is  taken  as  the  sum  of  the  three 
layer  thicknesses  d  =  t1  +  t2  +  tc»  assuming  that  the 
bonding  material  has  zero  thickness. 


The  faces  are  assumed  to  be  conventional  thin  plates 
with  orthotropic  material  properties  denoted  by  Eix»  Eiy»  Gixy' 
v.  ,  v.  ;  however,  computations  are  based  on  effective 


vixy'  viyx'  however'  computations  are  based  on  effective 

compressive  moduli  defined  by  ET  =  /z.  E"!  .  The  core  is 

x  3.x  x  y 

assumed  to  have  no  inplane  stiffness  at  all;  however,  the 


transverse  shear  moduli  are  represented  by  G  and  G  , 

cxz  cyz 

and  the  transverse  elastic  modulus  is  given  by  E 

c  z 


The  permissible  loadings  consist  of  uniform  inplane 
compression  loads  as  shown  in  Figure  3.  A  single  load 
condition  can  consist  of  a  uniform  load  N  compressing  the 
panel  in  the  x-direction,  or  a  uniform  load  N^  compressing 
the  panel  in  the  y-direction. 
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2.4  SANDWICH  PANEL  DESIGN  PARAMETERS 

For  the  purpose  of  this  sandwich  panel  optimization  study 
it  is  assumed  that  the  planform  dimensions  of  the  panel  and  the 
material  of  the  faces  and  core  are  predetermined.  The  design 
parameters  d1,...,dn  of  Equation  2.1  consist  of  the  thicknesses 
of  the  constituent  parts  of  the  sandwich  panel  together  with 
parameters  that  control  the  density  of  the  core  depending  on 
the  type  of  core  used  -  either  honeycomb  or  foam. 

For  the  design  of  sandwich  panels  having  honeycomb  core 
there  are  five  possible  design  parameters: 


1. 

Face  1  thickness, 

*1 

2. 

Face  2  thickness. 

t2 

3. 

Core  thickness, 

tc 

4. 

Honeycomb  cell  size, 

s 

5. 

Honeycomb  foil  thickness, 

In  the  case  of  foam  core  sandwich  four  design  parameters  are 
available: 


1. 

Face 

1  thickness. 

fcl 

2. 

Face 

2  thickness. 

t2 

3. 

Core 

thickness. 

fcc 

4. 

Core 

density. 

Pr 

The  various  possible  design  parameters  noted  above  can  be 
selected  in  any  numbers  and  combinations. 

2.5  SANDWICH  PANEL  MERIT  FUNCTIONS 

The  merit  or  objective  function  W(d^,...,dn)  of  Equation 
2.1  has  been  taken  to  be  the  total  weight  of  a  sandwich  panel. 

The  bonding  material  weight  will  be  approximately  constant  for 
a  panel  of  specified  planform  dimensions  regardless  of  the 
particular  values  of  the  design  parameters;  therefore,  the  weight 
of  the  bonding  agent  has  been  excluded  from  the  panel  weight  W. 
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The  expressions  for  the  objective  function  for  particular 
types  of  core  materials  are: 

1.  Hexagonal  Cell  Honeycomb  Core 

8tf  Pf 

W(t1,t2,tc,s,tf)  «  (pjtj^  +  P2t2  +  3^-  ■  tq)A  (2.3a) 

2.  Square  Cell  Honeycomb  Core 

2tf  Pf 

W(t1,t2,tc,s,tf)  =  (p^  +  p2t2  +  — g-^JA  (2.3b) 

3.  Foam  Core 

W(t1,t2,tc,pc)  =  (p^  +  p2t2  +  p  tfi)A  (2.3c) 

The  design  parameters  have  been  defined  in  the  previous  paragraph. 
The  face  and  honeycomb  material  densities  p^,  p2 ,  and  p^  are 
considered  to  be  predetermined  quantities,  as  is  the  panel 
pi an form  area  A. 


2.6  SANDWICH  PANEL  DESIGN  PARAMETER  CONSTRAINTS 

To  make  the  minimum  weight  design  of  sandwich  panels  more 
realistic,  numerous  design  parameter  constraints  are  imposed. 

All  of  the  constraints  are  of  the  inequality  type  (N  >  0,  M  =  0 
in  Equation  2.1). 

The  particular  inequality  constraints  g.  >_  0,  i  =*  1,...,N 
imposed  on  the  design  parameters  of  Paragraph  2.4  are  given  below. 
Each  constraint  is  designed  so  that  if  it  is  close  to  being 
violated,  then  the  constraint  value  will  be  a  small  positive 
number. 

1.  Minimum  Face  1  Thickness 


gx  = 


timin 


1 


(2.4) 
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2.  Minimum  Face  2  Thickness 


g2  =  t 


-1 


2  min 


3.  Maximum  Sandwich  Thickness 


g3  =  l  - 


*i  +  fc2  *  fec 
dmax 


4.  Minimum  Core  Thickness 


g4  =  t~ 
cmin 


-1 


5.  Minimum  Honeycomb  Foil  Thickness 


9c:  * 


5  tfmin 


-1 


6.  Minimum  Honeycomb  Cell  Size 


9*  * 


’6  s  ■ 
mxn 


-1 


7.  Minimum  Foam  Core  Density 

P„ 


g7  p 


-1 


cmin 


8.  Maximum  Honeycomb  Cell  Size 


max 


s8  s 


-1 


9.  Compressive  Yield  of  Face  1 


*9  *  1  ' 


N  E'/(tjE'  +  t2E') 
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(2.5) 


(2.6) 


(2.7) 


(2.8) 


(2.9) 


(2.10) 


(2.11) 


(2.12) 


10.  Compressive  Yield  of  Face  2 


g10 


=  1  - 


N  +  t2E') 


11.  Wrinkling  of  Face  1  -  Honeycomb  Core 

N  E{/{t1E£  +  ^E') 


gn  =  1  - 


(Flw> 


he 


where, 


(Flw>  =  (Equation  3:6,  Reference  6) 
he 


12.  Wrinkling  of  Face  2  -  Honeycomb  Core 

N  E'/ftjE^  +  t2E2) 


g12  *  1 


he 


where. 


(F2w)  =  (Equation  3:6,  Reference  6) 

he 


13.  Wrinkling  of  Face  1  -  Foam  Core 


*13 


N  EjVUjE'  +  t2E2 
—  1  —  - 


(Flw> 


foam 


where , 


(F1m)  =  (Equation  3:1,  Reference  6) 

AW  foam 


14.  Wrinkling  of  Face  2  -  Foam  Core 

N  E2/(tlE£  +  t2E2) 


*14 


=  1  - 


^7 


foam 
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(2.13) 


(2.14) 


(2.15) 


(2.16) 


(2.17) 


> 


Not  all  of  the  constraints  will  necessarily  be  used  for  a 
particular  design  problem.  For  example,  if  the  thickness  of 
Face  1  is  not  taken  as  a  design  variable  then  g^  would  not  be 
used.  In  the  case  of  multiple  loading  conditions,  multiple 
sets  of  constraints  9-18  are  necessary. 

The  constraint  functions  have  been  designed  so  that: 

(a)  a  positive  value  indicates  a  design  point  in  the 
acceptable  portion  of  the  design  space, 

(b)  a  negative  value  indicates  a  design  point  in  the 
unacceptable  region ,  and 

(c)  a  positive  value  near  zero  indicates  a  design 
point  on  the  boundary  of  the  acceptable  region. 
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SECTION  3 
IMPLEMENTATION 


This  section  contains  information  concerning  the  computer 
program  which  implements  the  procedure,  presented  in  Section  2,  for 
designing  sandwich  panels  for  minimum  weight.  The  computer 
program  SANOPT  is  operational  on  the  ASD  CYBER  computers  at 
Wright-Patterson  Air  Force  Base.  The  FORTRAN  source  code 
(Appendix  B)  is  resident  on  the  ASD  permanent  disc  storage 
files.  The  present  version  of  the  program  is  intended  for 
interactive  usage  through  remote  terminals  via  INTERCOM.  The 
program  operates  interactively,  and  executes  entirely  in-core. 
21,320  octal  words  of  central  memory  are  required  to  load,  and 
37,700  octal  words  are  necessary  for  execution. 


3.1  PROGRAM  ORGANIZATION 

SANOPT  is  organized  in  overlay  format  with  one  main 
overlay  and  five  primary  overlays.  The  overlay  structure  is 
as  follows: 


OVERLAY  OVERLAY 
LEVEL  NAME 


(0.0)  SANOPT 


(0.1)  LONGER 


(2.0)  SHORT 


(3.0)  SUMT 


OVERLAY 

FUNCTION 


Main  overlay  and  program  driver. 
Contains  the  field  interactive  input 
data  reader.  Reads  in  control 
information,  data  to  start  the 
problem  and  directs  calls  to  the 
primary  overlays. 

Primary  overlay  to  read  input  data 
and  initialize  design  parameters. 

The  parameters  along  with  their 
description  are  printed  out  on  the 
user  terminal. 

Primary  overlay  same  as  above,  but 
only  the  parameter  names  without  their 
description  are  printed  out  on  the 
user  terminal. 

Primary  overlay,  contains  the  program 
for  sequential  unconstrained 
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minimization  algorithm  (see  Reference  7 
for  details) . 

(4.0)  OUTPUT  Primary  overlay  contains  statements 

used  to  print  out  information  on  the 
results  of  the  optimization  problem. 

(5.0)  ALTER  Primary  overlay  contains  input  data 

reader.  It  allows  the  user  to  change 
any  data  read  in  the  overlays  LONGER 
or  SHORT. 


3.2  INPUT  VARIABLES 

The  following  describes  the  quantities  which  are  requested 
interactively  by  SANOPT: 


VARIABLE 

DEFINITION 

FCT 

Type 

of  Core 

=  Y: 

Honeycomb 

=  N: 

Foam 

TPE 

Type 

of  Cell  for  Honeycomb  Core 

=  H: 

Hexagon 

=  S: 

Square 

FI 

Type 

of  Face  Panels 

=  Y: 

Isotropic 

=  N: 

Orthotropic 

LTP 

Load 

Type 

=  1: 

Edgewise  compression  in  X-direction 

=  2:  Edgewise  compression  in  Y-direction 
=  3:  Edgewise  shear 
=  4 :  Edge  moment  in  X-direction 
=  5:  Edge  moment  in  Y-direction 
=  6:  Transverse  shear  in  X-direction 
=  7:  Transverse  shear  in  Y-direction 
=  8 :  Transverse  normal  load 
=  9:  Combined  loads 
=10:  Multiple  loads 

(Types  3-10  are  not  active) 

I ED  Panel  Edge  Conditions 

=  1:  All  edges  simply  supported 
=  2:  X-direction  clamped 
=  3:  Y-direction  clamped 
=  4:  All  edges  clamped 
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AA 

BB 

XNBAR 

NLT 

LC 

N 

T1,T2 

TC 

RHOC 

S 

TF 

IDV 

D 

RHO , RHOF 
RHOC 
E 
F 

GO 

R 

Ej  ,Ejx,Ejy 

PRiX,PRjY 

Gyxj 

RHOj 

YIELD j 

DFj 


Edge  Dimension  in  X-Direction 
Edge  Dimension  in  Y-Direction 
Applied  load 

Number  of  different  load  types 
Load  type 

Number  of  design  variables 
Thickness  of  face  panels 
Thickness  of  core 
Density  of  foam 

Size  of  cell  for  honeycomb  core 

Thickness  of  foil  for  honeycomb  core 

Design  variable  code 

Thickness  of  sandwich  panel 

Density  of  foil  for  honeycomb  core 

Density  of  foam 

Young's  Modulus  of  core 

Compressive  strength  of  core 

Transverse  shear  modulus  of  core 

Ratio  G  and  G^„  of  the  core 
cx  cy 

Young's  moduli  for  face  j 
Poisson's  ratio  for  face  j 
Shear  modulus  for  face  j 
Density  of  face  j 
Compressive  yield  stress 
Deflection  waviness  of  face  j 
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3.3  SAMPLE  TERMINAL  SESSION 

The  input  variables  defined  in  Paragraph  3.2  are  input 
interactively  by  the  user  of  SANOPT  in  response  to  prompts  by 
the  program.  The  format  of  the  prompts,  together  with  the 
variable  definitions  in  the  previous  paragraph,  make  the  data 
input  largely  self-explanatory.  This  paragraph  contains  a 
discussion  of  a  sample  terminal  session  which  illustrates  the 
use  of  SANOPT  to  determine  the  particular  configuration  of  a 
sandwich  panel  which  has  minimum  weight. 

The  objective  of  the  design  session  is  to  select  the 
thicknesses  of  the  two  faces,  the  thickness  of  the  hexagonal 
honeycomb  core,  the  honeycomb  cell  size,  and  the  honeycomb 
foil  thickness  of  a  sandwich  panel  with  characteristics  as 
shown  in  Figure  4.  The  simply  supported  panel  has  planform 
dimensions  of  80  inches  by  40  inches  and  is  subjected  to  an 
in-plane  compressive  load  of  1000  lb/ in  on  the  40  inch  side. 

The  upper  and  lower  face  sheets  are  aluminum  with  Young's 
modulus  30  x  106  psi ,  Poisson's  ratio  .25,  and  weight  density 
.1  lb/in3.  The  core  is  to  be  an  aerospace  grade  5052  alloy 
hexagonal  aluminum  honeycomb  selected  from  Figure  5. 

Figures  6a-l  contain  a  sample  interactive  session  with  the 
SANOPT  sandwich  panel  optimization  program.  In  the  '*&?v’le 
session  the  printouts  which  have  not  been  underlined  are  either 
informative  comments  or  prompts  which  require  some  action  by 
the  user:  the  underlined  printouts  are  the  responses  by  the  user 
to  the  various  prompts.  Each  of  the  Figures  6a-l  are  considered 
in  turn  below: 

(a)  Figure  6a 

The  "LOGIN"  procedure  and  the  program  access  routine 
are  illustrated  in  Figure  6a.  After  the  computer  (the  ASD 
CDC  system)  responds  to  the  telephone  dialup  by  identifying  the 
system,  the  date,  and  the  time,  the  user  responds  to  a  series 
of  prompts  by  typing  "LOGIN",  a  valid  problem  number,  a  password, 
and  a  terminal  identification  number.  The  login  process  is  then 
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5052  ALLOY  HEXAGONAL  ALUMINUM  HONEYCOMB 
AEROSPACE  GRADE 


H  EXCEL 
HONEYCOMB 
DESIGNATION 

Cell-Material-Gage 


COMPRESSIVE 

0M 

Stabilise* 

Strength 

psi 

Strength 

P*i 

Modului 

ksJ 

*YP 

min 

typical 

290 

21S 

75 

S4S 

405 

150 

910 

680 

240 

1470 

1100 

350 

2325P 

— 

900P 

215 

160 

55 

410 

300 

110 

720 

535 

195 

1130 

800 

285 

1600 

1180 

370 

135 

too 

34 

290 

21S 

75 

S2S 

90S 

145 

BIO 

400 

220 

1  130 

800 

285 

1470 

1 100 

350 

95 

70 

20 

17S 

130 

45 

340 

250 

90 

SOS 

370 

140 

690 

510 

190 

880 

660 

235 

1420 

1050 

340 

45 

20 

10 

95 

70 

20 

175 

130 

45 

270 

200 

70 

390 

285 

105 

4S5 

355 

135 

745 

535 

200 

1020 

750 

265 

KATE  SHEAR 


"L"  Direction 


Strength 

P*i 


"W"  Direction 


Strength 

psi 


min  typical 

90  22.0 


Figure  5.  Sample  Honeycomb  Data. 
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ASD  COMPUTER  CENTER  INTERCOM  5.0 
SYSTEM  CSA 

DATE  12/05/80  TIME  12.01.05. 

PLEASE  LOGIN 
LOGIN 

ENTER  PROBLEM  NUMBER- _ 

ggBggjgiii  ENTER  PASSWORD- 
ENTER  3— DIGIT  TERMINAL  ID- _ 


12/05/80 


LOGGED  IN  AT  12.01.46. 
WITH  USER- ID  WH 
EQUIP/PORT  16/017 


LOGIN 


CREATED  12/05/80  TODAY  IS  12/05/80 


COMMAND-  ATTACH,  P ,  SANQFJ ,  CY^ID^SONI ,  SN=AFFDL 
COMMAND-  FTN. I ~  P , L—PR , B~L  GO 

7.698  CP  SECONDS  COMPILATION  TIME 
COMMAND-  LGO 


FOLLOWING  IS  THE  MOST  RECENT 
VERSION  OF  T;  1C  SANDWICH  PANEL 
OPTIMIZATION  PROGRAM  USED  TO  COMPUTE 
THE  MINIMUM  WEIGHT  OP  A  SANDWICH 
PANEL  REQUIRING  SEVERAL  DESIGN 
VARIADLC3. 

AS  Or  IMCW  THE  PROGRAM  IS  SET  UP 
I'G  OPTIMIZE  SANDWICH  PANELS  WITH 
EITHER  A  P GAM  OR  HONEYCOMB  CORE  WITH 
AN  EDGEWISE  COMPRESSIVE  LOAD  IN  THE 
X  - D I REC r I ON  OR  Y - D I RECT I ON . 

THE  PROGRAM  WILL  BE  MODIFIED  TO 
ACCOMODATE  OTHER  CONDITIONS. 

(  JULY  15,  1980)  U.D.R.I. 


Figure  6a.  LOGIN  and  Program  Entry. 
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;  "  ’  ■  *  ■  ^  1  f  i  i\  1_>  W 1 i  i  i  dUtlL  OF  T  I M 1 h  <  I  GN  PI  vUL/Fv  AM 

i_  ThL  .  ur.M  UF  '.ME  1  NF'Ul  DEC 1  RED  ( Y  -YES » N -NO )  V'  -  -- .  - 

i  :•  >i  i  t'juE  i  i.ijI  it  •- uRE  Dt.  _•  i RED  ( V  - YEc » N— NO ) '?  -  *  •  -  —  : 

!•;.  h  HEXAGON  OR  SQUARE  CELL  DESIRED (H, S)?  - -  — ii 

AKL  ISOTROPIC  FACES  DES I  RED  <  Y- YES » N~NO )  ?  - s  _Y 

.  k  **  *■«  **#***-<  b  *  a  LOAD  AND  EDGE  CONDITIONS  **************!*** 

ENTER  LOAD  Tv PC  DESIRED 

1-  EDGEWISE  COMPRESSION  X-DIR 

2 - EDGEWISE  COMPRESSION  Y-DIR 

NUMBER . .  -  -  -  - - - -  -  - - :  1 

ENTER  TTiC  PANEL  EDGE  CONDITIONS  — 

1-ALL  EDGES  SIMPLE  SUPPORT 
2  X-DIR  CLAMPED,  Y-DIR  SIMPLE  SUPPORT 
S- X-DIR  SIMPLE  SUPPORT,  Y-DiR  CLAMPED 
4 -ALL.  EDGES  CLAMPED 

DOER  I  iil_  EDGL  DIMENSION  IN  THE  X-DIR  - - - s  GO 

uNTER  THE  EDGE  DIMENSION  IN  THE  Y-DIR  - -  - - - :  40 

LN1  £E\  'HE  Ai  T'D.  i  ED  LOAD  (LB/  IN)  - - - - — - :  1000 

&..***■»**»**»**  DESIGN  VARIABLES  ********************* 

THL  DESIGN  VARIABLES  FOR  THE  PANEL  ARE  AS  FOLLOWS: 

i.  --  T 1  2  -  T 2  3-TC  4-S  5-TF 

ENTER  THE  NUMBER  OF  DESIGN  VARIABLES  AND 

DESIGN  VARIABLE  CODES  IN  ASCENDING  ORDER  - :  512345 

ENTER  THE  INITIAL  VALUES  OF  THE  DESIGN  VARIABLE'S 

i-Tl  2 - T  2  3-TC  4-S  5-TF 


~  —  -  .  — *  —  —  —  — :  p  01  ■  0 1  ■  7  Si  .  l S  •  00 1 

ENTER  THE  MINIMUM  VALUE  FOR  T1  - TTT -.mi... 1  T005“ 

ENTER  THE  MINIMUM  VALUE  FOR  T2  - :  77!^ 

ENTER  T.tE  MINIMUM  VALUE  FOR  TO  - :  71 

ENTER  THE  MINIMUM  AND  MAXIMUM  VALUE  FOR  £  - :  .  125  . 375 

ENTER  THE  MINIMUM  VALUE  FOR  TF  - :  .0007 

ENTER  THE  MAXIMUM  THICKNESS  OF  THE  SANDWICH  PANEL - :  773 


*«•***•**  ***********  ***  CORE  PROPERTIES  ****************** ** * 


ENTER  THE  DENSITY  OF  THE  FOIL  - : 

ENTER  THE  MODULUS  OF  THE  CORE  - : 

ENTER  THE  COMPRESSIVE  STRENGTH  OF  THE  CORE  - - 

ENTER  THE  TRANSVERSE  SHEAR  MODULUS  OF  THE  CORE  - : 

ENTER  THE  RATIO  GCX/GCY  OF  THE  CORE  - ; 

********** **********  FACE  1  PROPERTIES  ******************** 

ENTER  YOUNG  'S  MODULUS  OF  THE  FACE - : 

ENTER  POISSON  S  RATIO  OF  THE  FACE - : 

ENTER  THE  DENSITY  OF  THE  FACE - : 

ENTER  THE  COMPRESSIVE  YIELD  STRESS  OF  THE  FACE  - : 

ENTER  THE  DEFLECTION  WAV  I  NESS  OF  THE  FACE - : 

********************  FACE  2  PROPERTIES  ******************** 

ENTER  YOUNG'S  MODULUS  OF  THE  FACE  - : 

ENTER  POISSON'S  RATIO  OF  THE  FACE  - : 

ENTER  THE  DENSITY  OF  THE  FACE  - : 

ENTER  THE  COMPRESSIVE  YIELD  STRESS  OF  THE  FACE  - : 

ENTER  THE  DEFLECTION  WAV  I  NESS  OF  THE  FACE  - : 

ARE  THERE  ANY  CHANGES  TO  THE  DAT  A ( Y - YES . N-NO ) ? - : 


._1 

75E3 

290 


10E6. 


.  1 

50000 

.001 


10E6 


.  1 

goooo 


Diyure  6b.  Input  Data  for  Initial  Design  Iteration. 
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OPTIMIZATION  BEGINS 


PANEL  NO.  1 


******************  OPTIMIZED  SANDWICH  PANEL  -it*************** 

DESIGN  PARAMETERS 


THE  THICKNESS  OF  FACE  1  SHOULD  BE  - :  .024897 

THE  THICKNESS  OF  FACE  2  SHOULD  BE - :  .024897 

THE  THICKNESS  OF  THE  CORE  SHOULD  BE - :  1.042736 

THE  CELL  SIZE  OF  THE  CORE  SHOULD  BE - :  .245888 

THE  THICKNESS  OF  THE  FOIL  SHOULD  BE  - :  .004147 

THE  ABOVE  VALUES  OF  THE  DESIGN  PARAMETERS  YIELDS 

A  MINIMUM  WEIGHT  FOR  THE  SANDWICH  PANEL  TO  BE  - :  30.939994  LBS 


THE  FINAL  CONSTRAINT  VALUES  ARE: 

MINIMUM  FACE  1  THICKNESS  . 39794352E+01 

MINIMUM  F ACC  2  THICKNESS  . 39794326E+01 


MAXIMUM  SANDWICH  THICKNESS 
MINIMUM  CORE  THICKNESS 
COMPRESSIVE  YIELD  FACE  1 
COMPRESSIVE  YIELD  FACE  2 
MINIMUM  FOIL  THICKNESS 
MINIMUM  CELL  SIZE 


. 27 1 64630E+00 
. 9427362 1E+01 
. 59834792E+00 
. 59s S4792E+00 
-  A  V  >z»  1E+01 
. 9671 037 1E+00 


WRINKLING  FACE  I 
WRINKLING  FACE  2 
DIMPLING  TACE  1 
DIMPLING  FACE  2 

BUCKLING  LOAD 
MAXIMUM  CELL  SIZE 


. 74174157E-06 
. 21017850E-06 
. 908 1 S034E+00 
. 908 1 S024E+00 
. 1 1335903E-06 
. 32308482E+00 


Figure  6c.  Initial  frial  Design. 
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>N  CO  IP  -0 


ARE  THERE  ANY  CHANGES  TO  THE  DATA ( Y- YES, N-NO ) ?  - 

IS  A  LIST  OF  POSSIBLE  CHANGES  DESIRED < Y-YES, N-NO ) ?  - 

I- TYPE  OF  HONEYCOMB  CORE 
-EDGE  CONDITIONS 
-DIMENSION  X-DIR 
••DIMENSION  Y-DIR 
-APPLIED  LOAD 
-INITIAL  VALUES 
-MINIMUM  T 1 

8- MINIMUM  T"2 

9- MINIMUM  TC 

1 0-  M I N I MUM  RHOC ( FOAM ) 

I I - MINI MUM  S ( HONEYCOMB ) 

12- MAXIMUM  S ( HONEYCOMB ) 

1 3-  MINI MUM  TF <  HONEYCOMB ) 

14- T  1 

15-  T  2 

io-rc 

1  7  -RHOC ( FOAM ) 

18- 3 (HONEYCOMB) 

19-  TFT HONEYCOMB) 

20 - CORE  MATERIAL  DENSITY 
21  -CORE  MODULUS 

22-  CORE  COMPRESSIVE  STRENGTH 

23 - CORE  TRANSVERSE  SHEAR  MODULUS 

24-  CORE  OCX/GCY  RATIO 

25-  VOUNGS  MODULUS  FACE  1 
26  ■YOUNGS  MODULUS  FACE  2 
27 -POISSONS  RATIO  FACE  1 
2C-PGI3S0NS  RATIO  FACE  2 

29 - FACE  1  DENSITY 

30- FACE  2  DENSITY 

31 - FACE  1  YIELD  STRESS 
3i'->  ACE  L  YIELD  STRESS 
36  PACE  l  DEFLECTION 
34-FACE  2  DEFLECTION 

ENTER  THE  NUMBER  OF  CHANGES  AND  THE  CHANGE  CODES 

r  RUM  THE  LIST  - - - 8  6_ 

X(l--U)  . . —  .01  . 01  .  75  .  35  .  00 1 

ARE  THERE  ANY  OTHER  CHANGES < Y-YES, N-NG ) ?  - 


Figure  6d.  List  of  Possible  Changes  and  Data  for  New 
Starting  Point. 
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******* ************  OPTIMIZATION  BEGINS  ******************* 


FANEL  NO. 


*****************  OPTIMIZED  SANDWICH  FANEL  **************** 

DESIGN  PARAMETERS 


THE  THICKNESS  OF  FACE  1  SHOULD  BE - :  .024897 

THE  THICKNESS  OF  FACE  2  SHOULD  BE  - :  .024897 

THE  THICKNESS  OF  THE  CORE  SHOULD  BE  - :  1.042746 

THE  CELL  SIZE  OF  THE  CORE  SHOULD  BE  - .247551 

THE  THICKNESS  OF  THE  FOIL  SHOULD  BE  - :  .004175 

THE  ABOVE  VALUES  OF  THE  DESIGN  PARAMETERS  YIELDS 

A  MINIMUM  WEIGHT  FOR  THE  SANDWICH  PANEL  TO  BE  - :  30.939993  LBS 


THE  FINAL  CONSTRAINT  VALUES  ARE: 
MINIMUM  FACE  I  THICKNESS 
MINIMUM  FACE  2  t HICKNESS 
MAXIMUM  SANDWICH  THICKNESS 
MINIMUM  CORE  THICKNESS 
COMFRESSIVE  YIELD  FACE  1 
COMPRESSIVE  YIELD  FACE  2 
MINIMUM  FOIL  THICKNESS 
MINIMUM  CELL  SIZE 
WRINKLING  FACE  1 
WRINKLING  FACE  2 
DIMPLING  FACE  1 
DIMPLING  FACE  2 
DUCKLING  LOAD 
MAXIMUM  CELL  SIZE 


.  397947 28E+01 
. 39794740E+0 1 
.271 6398 1E+00 
. 94274554E+0 1 
. 598351 10E+00 
. 598351 10E+00 
. 49638305E+01 
. 98040829E+00 
. 23058593E— 06 
.48491 06 IE-06 
. 90693623E+00 
. 90693628E+00 
. 949 1 9045E-07 
.5148391 5E+00 


Figure  6e.  Initial  Trial  Design  Recomputed. 
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ARE  THERE  ANY  CHANGES  TO  THE  DATA  ( Y -YES.  N-NG )  ?  - :  _Y 

IS  A  LIST  OF  POSSIBLE  CHANGES  DESIRED! Y- YES. N-NO)?  - :  N 

ENTER  THE  NUMBER  OF  CHANGES  AND  THE  CHANGE  CODES 
FROM  THE  LIST  - :  1  6 

XC1-N)  - .01  .  01  .  25  .  25  .  00 1 

ARE  THERE  ANY  OTHER  CHANGES <Y- YES. N-NO)?  - s  N 


**#’H-*-i**'*#**#tt*#***it  OPTIMIZATION  BEGINS  ***■«#*****•**•*#■*#*** 


PANEL  NO. 


****■«•*•*#«•■««•»(••»**»*  OPTIMIZED  SANDWICH  PANEL  #********•#■***#*# 

DESIGN  PARAMETERS 


THE  THICKNESS  OF  FACE  1  SHOULD  BE  - :  .024900 

THE  THICKNESS  OF  FACE  2  SHOULD  BE  - s  .024900 

THE  THICKNESS  OF  THE  CORE  SHOULD  BE  - - -  1 . 042857 

THE  CELL  SIZE  OF  THE  CORE  SHOULD  BE  - s  .  240805 

THE  THICKNESS  OF  THE  FOIL  SHOULD  BE  - - .004146 

THE  ABOVE  VALUES  OF  THE  DESIGN  PARAMETERS  YIELDS 

A  MINI MUM  WEIGHT  FOR  THE  SANDWICH  PANEL  TO  DE  - :  30.939994  LBS 


FINAL  C 

ONSTRAINT 

VALUES  ARE: 

MINI MUM 

F ACC  1  TH 

1 CKNESS 

. 39799470E+0 1 

MINI MUM 

FACE  2  IH 

I  CKNESS 

. 39799502E+0 1 

MAXIMUM 

SANDWICH 

T H I CKNESS 

. 27 1 56223C+00 

MINI  i .  JM 

CORE  MIC 

KNESS 

. 942857 16E  <-01 

compress 

'i’v'E  V 1  Li_D 

FACE  1 

.  5 9 0 3 0 9 4 3 E + 0  0 

COMPRESS 

I VE  YIELD 

FACE  2 

.  5 9 8 3 8 9 4 3 E  +•  0  0 

1  •;  1  N  I  MUM 

FOIL  THUS 

KNESS 

. 49224409C +0 1 

MINI  MUM 

CELL  SIZE 

. 9670765 1E+00 

WRINKLING  1  ACL  1 

. 25608 122E-06 

WRIMKlING  FACE  2 

. 889 1 0409 E -06 

I.  1  T1PL I NG 

!  Act.  1 

. 90021 I24E+00 

DIMPLING 

FACE  2 

.  '90821 1  SOL +00 

DUCKLING 

i  LOAD 

.  4 8 9 0 8 5 5 3 E  -  0  7 

MAX  T  Ml  IM 

r-Fi  i  t  ? r 

. 52510590E+00 

Figure  6f.  Data  for  New  Starting  Point  and  Recomputed 
Initial  Trial  Design. 
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ARE  THERE  ANY  CHANGES  TO  THE  DATA < Y- YES, N-NG) ?  - : 

IS  A  LIST  OF  POSSIBLE  CHANGES  DESIRED < Y -YES, N-NO > ? - : 

ENTER  THE  NUMBER  Or  CHANGES  AND  THE  CHANGE  CODES 

FROM  THE  LIST  - :  4  21  22  S 

R  . . . . . 2.46 

ARE  THERE  ANY  OTHER  CHANGES ( Y -YES , N-NO ) ?  - : 

*•*«*•*•*•* » « « « fc »  OPTIMIZATION  BEGINS  »•»•#**•*#«•****# 


23  24 


****  h -rf**  *  a  *  OPTIMIZED  SANDWICH  PANEL  a  a  «  a  a  *  s  a  «  »*»***■*• 
DESIGN  PARAMETERS 


THICKNESS  OF  FACE  1  SHOULD  BE  - 

rKICKNLSS  OF  FACE  2  SHOULD  BE  — 

TK i  -KNESS  Or  •  hL  CORE  SHOULD  DE  -- 

CEl..i.  SIZE  OF  7::~  CORE  SHOULD  DE  — 

F  i-.  LSKNLSS  -■!  1 4  .i_  F C I L  SHOULD  DE  - — 


.018335 
. 018355 
.957180 
. 235426 
. 0020 1 1 


HIE  DESIGN  PARAMETERS  YIELDS 
’HE  SANDWICH  F ’ANEL  iU  DL  -  - 


18. 724312 


L1*  I  t".  ■  }  I  .-*•  »“  r-,  -  ,  4  ~r  I  «  ,"\  I  I  trr 

>  l.'ti-iL.  C'  -•!  To  i  iMiii'i  I  v  ni\u  • 

MI.  4 1  MUM  i  t.CC  1  .  i.lcKNESS 
M  IN.  MUM  ;  ,  .ME  2  K -SICKNESS 


i'  1 1—.  *  1 1  i  Ul  I  i.ANDW  i  i ;  T  HI  C 
T.  t  A  i  4>. IM  ■.  OR.-  rHiCKNES 
.  I ■  '•  I  72  FiLl-D  ;  AC 
Cui-,1  MESS  I  VC  YIL..D  FAC 
.  •  i  N I  ' .  Ji  I  f  .j  A  L  THJ NLS 
MINIMUM  CELL  SIZE 
WK  i.  i  J  i ,  C I  i'6  j  F  AC  E  1 

wrinkling  face  2 

D 1 1 .  i  Nu  ;  Act-  I 
D  1 I NU  I  AtT.  2 
B  1 1  <U  _CAD 

MAXIMUM  CELL  SIZE 


KNESS 


.  267 1 07 06E+0 1 
.  267 106.5  6E-+- 01 
.  S3740614C+00 
,  857 1 3009L  +  0 1 
,  4551 9973C+00 
,  4351 9 970 E + 0 0 
,  1 8728077E  +  0 1 
,  88341 166E  +  00 
,  7914971  IE- 06 
,51 28903 IE-06 
,  78994470E+00 
,  7 8 9 9 4  4 5 8 E + 0  O 
,  797891 96C -07 
,  5928541 1E+00 


Figure  6g.  Change  Data  for  1/4  -  5052  - 
Associated  Design  Results. 


004  Core  and 


ARE  THERE  ANY  CHANGES  TO  THE  DATA  <.  Y -  YES,  N  -  NO )  '•  - - :  _Y 

IS  A  LIST  or  POSSIBLE  CHANGES  DEC I RED ( Y  YES  >  N  -NO ) ?  -  - . :  N_ 

ENTER  THE  NUMBER  01  CHANGES  AND  THE  CHANGE  CODES 

FROM  HIE  LIST  . . ■  ■ . —  - .  -  4.  21  22  Si  Sj 

E  . . -  - 140ES 

GO  - - -  -6GE3 

ARE  THERE  ANY  OTHER  CHANGES  <  Y -YES » N-NO ) ?  - - :  N_ 

*•»•»*•»•» »•*•***•#■»•« ***#• *  OPTIMIZATION  BEGINS  ■***■**•#*•»*#•*#•**#**** 


PANEL  NO. 


*■»■#*■«■  tt -a  #•#■**■#•»  •«•*#*  OPTIMIZED  SANDWICH  P  ANEL  •»#»•#*«*#  *  ##*•»##  # 

DESIGN  PARAMETERS 


THE  THICKNESS  OF  FACE 
THE  THICKNESS  OF  FACE 


SHOULD  BE  - 

SHOULD  DE  - 


1  HE  CHIC 


KNES- 


THE  CELL  SIZE  OF 
THE  ! HICKNESS  OF 


THE  CORE 
THE  CORE 
THE  FOIL 


SHOULD  DE 
SHOULD  BE 
SHOULD  BE 


1.010501 
. 2 4 EC Sc 
. 00 S i 20 


THE  ABOVE  VALUES  OF  THE  DESIGN  PARAMETERS  YIELDS 
A  MINIMUM  WEIGHT  FOR  THE  SANDWICH  PANEL  TO  DE  -  - . :  25.374446 


THE  FINAL  CONSTRAINT  VALUES  ARE: 

MINIMUM  PACE  i  THICKNESS  . 34271 442E +01 

MINIMUM  FACE  2  THICKNESS  .34271 CO IE +01 


MAXI MUM  SANDN i Cl  I  TH I CKNESS 
MINIMUM  •-•URL  ■  H  IlKNEsS 
COMi'RESSIVE  YIELD  FACE  I 


.  291431 9 1C+00 
.91 858067 L~-  '  i 
. 54324 182E+0G 


COMPRESS 1 VC  YIELD  FACE  2  . 54024 1S2E+00 
M I N I MUM  FO i L  1 H I CKNESS  . 34685399E+0 1 
MINIMUM  CELl.  SIZE  .  9406S993E+00 


WRINKLING  FACE  I 
WRINKLING  FACE  2 
DIMPLING  FACE  i 


. 1 94  7 2559C-06 
. 155S7763E-05 
. S7283655E+00 


DIMPLING  r  AlE  2 
BUCKLING  LOAD 
MAXIMUM  CELL  SIZE 


. S 7 2 S 3 6 9 0 E+00 
. 2 1 0742 39C-06 
.  54584200E+00 


Figure  6h.  Change  Data  for  1/4  -  5052  -  .002  Core  and 
Associated  Design  Results. 
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ARE!  (HERE  ANY  CHANGES  TO  THE  DATA ( Y-YES.  N-NO )  ?  -  - :  .V 

IS  A  LIST  OF  POSSIBLE  CHANGES  DES I RED (Y-YES, N-NO >?  - :  N 

ENTER  (HE  NUMBER  OF  CHANGES  AND  THE  CHANGE  CODES 

FROM  THE  LIST  - - - - :  4  21  22  23  24 

E  - 90E3  - 

F  - - - - 340 

GO  — . — - 50E3 

R  - - --2.  OS 

ARE  THERE  ANY  OTHER  CHANGES < Y-YES , N-NO ) ?  - :  N 


*******************  OPTIMIZATION  BEGINS  ******************* 

PANEL  NO.  6 


*****************  OF'TIMIZEEi  SANDWICH  PANEL  ****** ******** ** 

DESIGN  PARAMETERS 


THE  THICKNESS  OF  FACE  1  SHOULD  BE  - :  .024058 

THE  THICKNESS  OF  FACE  2  SHOULD  BE  - :  .024058 

THE  THICKNESS  OF  THE  CORE  SHOULD  BE  - - - :  1.036955 

THE  CELL  SIZE  OF  THE  CORE  SHOULD  BE  - :  .249145 

THE  THICKNESS  OF  THE  FOIL  SHOULD  BE - :  .  003S90 

THE  ABOVE  VALUES  OF  THE  DESIGN  PARAMETERS  YIELDS 

A  MINIMUM  WEIGHT  FOR  THE  SANDWICH  PANEL  TO  BE  - :  29.213107  LBS 


THE  FINAL  CONSTRAINT  VALUES  ARE: 
MINIMUM  FACE  1  THICKNESS 
MINIMUM  FACE  2  THICKNESS 
MAXI  MUM  SANDW I CH  TH I CKNESS 
MINIMUM  CORE  THICKNESS 
COMPRESSIVE  YIELD  FACE  1 
COMPRESSIVE  YIELD  FACE  2 


.381 16619E+01 
. 33116575E+01 
.27661922E+00 
.  93695458E+0 1 
.  58434 300E+00 
.  58434300E+00 


MINIMUM  FOIL  THICKNESS 
MINIMUM  CELL  SIZE 
WRINKLING  FACE  1 
WRINKLING  FACE  2 
DIMPLING  FACE  1 
DIMPLING  FACE  2 
BUCKLING  LOAD 
MAXIMUM  CELL  SIZE 


. 45571513E+01 
. 993 1 6288E+00 
. 15994002E-05 
. 683277 17E-06 
. 395522 /4E+00 
. 89552255E+00 
. 33754923E-06 
. 5051 4543E+00 


Figure  6i.  Change  Data  for  1/4  -  5052  -  .0015  Core 
and  Associated  Design  Results. 
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ARE  THERE  ANY  CHANGES  TO  THE  DATA <  Y- YES >  N-NO ) ? - : 

IS  A  LIST  or  POSSIBLE  CHANGES  DESIRED < Y- YES, N- NO > ?  - : 

ENVER  THE  NUMBER  OF  CHANGES  AND  THE  CHANGE  CODES 
FROM  THE  LIST  - :  4_  21  22  23  24 

E  . . . . . -19oE3 

GO  — . . S2E3 

R  .  -—2.31 

ARE  T HERE  ANY  OTHER  CHANGES ( Y- YES. N-NO)?  - - 

« uit «  in*  **■*•«■  a-aa-a  -tt  *#-aa  OPTIMIZATION  BEGINS  -itaa-a-aa  a*  a  a  a  a  aa  a  a#*« 


PANEL  NO. 


*#«#■»■ a-  *  «  ««#««  *•»  OPTIMIZED  SANDWICH  PANEL  #*aa;Htaaaa«aas**- 

DESIGN  PARAMETERS 

THE  THICKNESS  ■  OF  FACE  1  SHOULD  BE  - - .020331 

THE  THI CKNLSS  OF  F  ACL  2  SHOULD  BE  - —  •  .  020SS1 

THE  VHICKNESS  OF  THE  CORE  SHOULD  BE  - - - :  .999022 

THE  CELL  SIZE  OF  THE  CORE  SHOULD  BE  - - .240961 

THE  THICKNESS  OF  THE  FOIL  SHOULD  BE  — - :  .002700 


1  HE  ABOVE  VALUES  OF  THE  DESIGN  PARAMETERS  YIELDS 

A  MINIMUM  WEIGHT  FUR  THE  SANDWICH  PANEL  TO  BE  - :  22.900492  LBS 


THE  FINAL  CONSTRAINT  VALUES  ARE: 

MINIMUM  FACE  1  THICKNESS  . 31661007E+01 


M I N I  MUM 

‘ACL  2  THICKNESS 

.  SI 661 01  4E-*  01 

M  A  X  j.  i  i'M  i 

..ANDW I CH  TH 1 CKNESS 

. 3062 1 100E+00 

MINI  MUM 

sure  i  h  i  lt.ness 

. 89902 24SEM 0 1 

COMPRESS 

IVE  YIELD  FACE  1 

.  5 1 9 9 3 4 8 3 E + 0  0 

CO  Ml  'REGS 

1VC  YIELD  FACE  2 

. S1993433C+00 

M I N l MUM 

OIL  TH I Cl, NESS 

. 236334  G5E  +  0 1 

i't  1 N 1  MUM 

ZELL  SIZE 

. 9276S639E+00 

WR  I  NFL  j.  N 

S  FACE  1 

. 36941077E  -06 

WRING.  ING  FACE  2 

. 53460600E-06 

D 1  MV  L.  j.  NO 

FACE  1 

. 34944 169E+00 

DiMPuiNG 

FACE  2 

. 84944 174E+ 00 

BUG  H  ING 

LOAD 

.  1 60763 2 4 E.  -06 

MAXIMUM 

CELL  SIZE 

. 55626974E+00 

Figure  6j.  Change  Data  for  1/4  -  5052  -  .0025  Core  and 
Associated  Design  Results. 
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>|zl 


690  S2E3  2. 31 
IF  - .  0025 

£  1  .  £2 »  PR  I  .  PR,. ,  RHO 1 ,  RH02 ,  DF 1 ,  DF2  —  1 0E6  1 0E6  .25  .  25  ._1  ^ 
COMPRESSIVE  STRESS  VIELD1 , YIELDS  — 50000  50000 

ARE  THERE  ANY  CHANGES  TO  THE  DATA ( Y -YES , N-NQ ) ?  . . . 

** » ii  ft  shh*  •**  -ihi-  *  OPT  I  MI  ZATION  BEGINS  **■  »##*»»#*#»*  ft  -s  mum 


PANEL  MO. 


i 


ft-mm-ft  ft-s  mum  -ft  -ft-ft -mt  ft  uF’TIMI  ZED  SANE'UIICH  PANEL  #-ft  im-im  sm-M  im  ft  ft  fttm 

DESIGN  PARAMETERS 


THE  THICKNESS  OF  FACE  1  SHOULD  DE  - 

THE  THICKNESS  OF  FACE  2  SHOULD  DE  - 

THE  THICKNESS  UF  THE  CORE  SHOULD  BE  --- 


. 021904 
•  O  i  0  4 
1 . 050  i  9  0 


THE  ABOVE  VALUES  OF  THE  DESIGN  PARAMETERS  YIELDS 
A  MINIMUM  WEIGHT  FOR  THE  SANDWICH  PANEL  TO  DE  - : 


uBS 


THE  FINAL  CONSTRAINT  VALUES  ARE: 
MINIMUM  FACE  i  THICKNESS 
MINIMUM  FACE  2  THICKNESS 
MA X  X I  iUiv1  SANDW I  CM  TH I CKNESS 
M 1 N i MUM  CORE  TH I CKNESS 
COMPRESSIVE  YIELD  FACE  1 
COMPRESSIVE  YIELD  PACE  2 

wrinkling  pace  ; 

WRINKLING  PACE  2 
Ell  Ml  ’LING  rACE  1 

DIMPLING  FACE  2 
BUCKLING  LOAD 


.  33008869E  KU. 
. 33308077E+0 1 
. 27066224E+0C 
. 9501 9776E+01 
. 54347 14SE+00 
. 54347 14GE+00 
. 296666 14E-06 
. 4639369 IE -06 
. 36062 17 1 Ei 00 
. 36062 I76E+C0 
.  82491 468E -07 


Figure  6k.  Short  Form  of  Data  Input,  Fixing  the  Core 
Properties,  and  Final  Design  Results. 
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IzHIz 


ARE  THERE  ANY  CHANGES  TO  THE  DATA< Y-YES, N-NO) 

IS  A  NEW  PROBLEM  DESIRED (Y-YES, N-NO)?  - 

STOP 

037700  MAXIMUM  EXECUTION  FL. 

12.039  CP  SECONDS  EXECUTION  TIME. 
COMMAND-  LOGOUT 

CPA  26.  £53" "SEC.  16.529  ADJ. 

10  54.409  SEC.  16.104  ADJ. 

CRUS  39.359 

CONNECT  TIME  0  HRS.  19  MIN. 

11/21/30  LOGGED  OUT  AT  12.19.55. 


Figure  6k.  (concluded) 


complete.  A  series  of  commands  is  then  issued  to  retrieve  the 
FORTRAN  source  code  from  permanent  file  storage,  to  compile  the 
source  code  and  create  a  binary  deck,  and  to  execute  the 
created  binary  code.  A  printed  message  indicates  that  the 
SANOPT  program  has  been  entered  successfully,  and  an  interactive 
design  session  has  been  initiated. 

(b)  Figure  6b 

Figure  6b  contains  a  series  of  program  prompts  and 
user  responses  to  define  completely  an  initial  panel 
optimization  problem.  The  data  are  input  in  seven  distinct 
sets  as  designated  below: 

1 .  General  Data 

The  first  group  of  input  data  prompts  and 
responses  is  a  series  of  questions  which  have  non-numerical 
responses.  There  are  two  forms  of  data  input,  the  long  form 
and  the  short  form.  In  the  case  of  the  long  form,  the  prompts 
are  self-explanatory.  The  short  form  generally  takes  less 
time  for  data  input,  but  the  prompts  are  brief,  with  the  data 
required  being  identified  only  by  the  variable  names  of 
Paragraph  3.2.  The  long  form  is  requested  here;  the  short 
form  will  be  considered  later.  A  honeycomb  type  core  with 
hexagonal  cells  is  specified  and  the  faces  are  identified 
to  be  isotropic. 

2 .  Loading  and  Edge  Conditions 

The  load  type  is  specified  to  be  an  edgewise 
compressive  load  on  the  x-edges  of  the  panel,  the  panel  is 
simply-supported  all  around,  the  planform  dimensions  are 
80"  x  40",  and  the  applied  edgewise  compression  load  is 
1000  lb/in. 

3.  Design  Variable  Data 

The  possible  sandwich  panel  design  variables 
are  defined  in  Paragraph  2.4.  In  the  case  of  honeycomb  core 
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these  are:  the  thickness  of  face  1  (Tl),  the  thickness  of  face  2 
(T2),  the  depth  of  the  core  (TC) ,  the  honeycomb  cell  size  (S) , 
and  the  honeycomb  foil  thickness  (TF) .  The  user  is  free  to 
define  any  or  all  of  these  possible  design  variables  as  the 
actual  design  variables.  In  this  case  all  five  of  the 
possibilities  are  selected  as  design  variables.  The  initial 
values  of  the  design  variables  define  a  point  in  the 
multi-dimensional  design  space  from  which  the  optimization 
process  will  start.  The  parameters  in  those  constraints  which 
apply  directly  to  the  design  variables  (see  Paragraph  2.6)  are 
then  defined  by  specifying  minimum  and/or  maximum  values  of  the 
design  variables.  Here  the  minimum  gage  of  the  faces  is  taken 
to  be  .005"  the  minimum  core  thickness  is  .1",  the  minimum  and 
maximum  cell  sizes  are  1/8"  and  3/8",  the  minimum  foil  gage  is 
.0007",  and  the  maximum  total  sandwich  is  defined  to  be  1.5". 

4 .  Core  Data 

The  next  set  of  input  data  refer  to  the  sandwich 
core.  In  our  example,  the  core  is  aluminum  with  weight  density 
of  .1  lb/in^.  The  transverse  modulus,  compressive  strength, 
shear  modulus  (in  the  load  direction),  and  the  ratio  of  the  two 
shear  moduli  have  been  taken  arbitrarily  (for  the  time  being) 
to  correspond  to  the  honeycomb  core  1/8  -  5052  -  .0007  in 
Figure  5.  These  values  must  be  set  before  the  minimization 
is  initiated  because  they  are  used  in  the  computation  of  some 
of  the  constraints  (Paragraph  2.6).  Obviously,  the  minimum 
weight  core  configuration  will  not  in  general  correspond  to  the 
assumed  initial  core  properties.  Our  procedure  will  be  to 
determine  a  minimum  weight  panel  using  a  core  with  properties 
in  Figure  5.  This  will  be  an  iteration  process  as  seen  below. 

5 .  Face  1  Data 

The  next  set  of  data  in  Figure  6b  refers  to 

face  1  properties.  In  this  case,  face  1  is  an  aluminum  sheet 

with  modulus  of  10  x  10^  psi,  Poisson's  Ratio  of  .25,  weight 

3 

density  of  .1  lb/in  ,  and  compressive  yield  stress  of  50,000  psi. 
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The  initial  waviness  of  the  face  sheets  (i.e.,  the  face  is  not 
perfectly  flat)  is  .001";  this  value  is  used  to  calculate  the 
panel  wrinkling  stress  (Chapter  3  of  Reference  6)  for  evaluating 
the  wrinkling  constraints. 

6.  Face  2  Data 

This  data  is  the  same  as  that  for  face  1  since  we 
have  identical  faces. 

7.  Data  Termination 

The  final  prompt  is  an  inquiry  as  to  whether 
there  are  any  changes  to  be  made  to  the  data.  At  this  point,  the 
answer  would  be  "Y"  if  some  mistake  had  been  made  in  entering 
any  of  the  data  in  Figure  6b;  then  the  user  would  be  permitted 
to  change  any  of  the  data.  In  this  case  we  choose  to  accept  the 
input  data  by  answering  "N".  This  causes  the  SANOPT  program 
to  terminate  the  data  input  phase  and  to  begin  the  optimization 
process. 

(c)  Figure  6c 

This  figure  presents  the  results  of  the  panel 
optimization  using  the  physical  data  and  the  initial  panel 
design  variables  provided  in  Figure  6b.  The  final  values  of 
the  design  parameters  are  indicated  as  well  as  the  minimum 
weight  computed  from  the  final  design  parameters.  Also  given 
are  the  values  of  the  various  constraint  functions  (Paragraph 
2.6)  which  apply  to  this  design  example.  In  this  case,  the 
final  values  of  the  design  parameters  have  been  influenced  by  the 
wrinkling  constraints  and  the  panel  buckling  constraint.  That 
is,  if  a  panel  were  built  with  the  calculated  design  parameters, 
it  would  be  that  acceptable  panel  with  the  least  weight  but 
would  be  on  the  verge  of  face  wrinkling  and  gross  panel  buckling. 

(d)  Figure  6d 

This  figure  indicates  the  procedure  for  making  one 
or  more  changes  to  the  set  of  data  input  in  Figure  6b.  Here  we 
specify  that  we  wish  to  change  data,  and  requea.  that  a  listing 
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of  all  possible  changes  be  presented;  in  response  a  list  of  34 
possible  data  changes  is  printed.  We  specify  that  a  single  data 
change  be  made,  and  in  particular  that  the  initial  values  of  the 
design  parameters  be  altered.  The  particular  values  of  the 
initial  parameter  values  are  then  input  for  the  five  design 
variables.  An  indication  that  no  other  data  changes  are  to 
be  made  then  initiates  the  optimization  process  with  the  new 
starting  point  in  the  design  space. 

The  reason  for  making  the  change  in  the  starting 
point  in  the  design  space  at  this  stage  of  the  panel  design 
process  is  to  ensure  that  a  valid  minimum  weight  has  been 
obtained,  and  not  just  a  relative  minimum.  Therefore,  a  new 
starting  point  is  specified  and  the  minimum  weight  is 
recomputed;  if  the  same  final  values  of  the  design  parameters 
and  the  same  minimum  weight  are  obtained  then  one  can  be 
relatively  sure  that  the  proper  minimum  has  been  found. 

(e)  Figure  6e 

This  figure  presents  the  results  of  the  optimization 
using  the  new  starting  point  in  the  design  space.  The  results 
are  the  same  as  those  computed  before. 

(f )  Figure  6f 

This  figure  defines  yet  another  design  starting 
point  but  this  time  not  requesting  a  complete  listing  of  the 
possible  changes.  The  results  again  are  identical  to  the 
original  problem. 

(g)  Figure  6g 

Although  the  design  computed  above  represents  a  valid 
optimum  design  for  the  data  input  in  Figure  6b,  the  core  is  not 
realistic  according  to  what  is  available  commercially.  The 
design  on  Figure  6c  calls  for  a  core  with  cell  size  of  about  1/4" 
and  a  foil  thickness  of  about  .004".  It  is  apparent  from  Figure 
5  that  no  core  is  available  with  these  dimensions  which  also  has 
the  material  properties  input  in  Figure  6b,  Therefore,  we  must 
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try  various  cores  from  Figure  5  until  we  get  a  minimum  design 
with  consistent  material  properties  and  geometry. 


Figure  6g  contains  alterations  to  the  data  to 
correspond  to  core  1/4  -  5052  -  004  in  Figure  5.  This  time  the 
optimized  design  calls  for  a  core  with  cell  size  again  1/4" 
but  with  foil  thickness  of  .002".  The  optimized  panel  weight 
is  some  12  lb.  lighter  than  the  weight  of  the  initial  trial 
design  but  it  is  not  feasible  because  there  is  no  core  available 
with  the  calculated  cell  size  and  foil  thickness  which  has  the 
properties  specified  in  Figure  6g. 

(h)  Figure  6h 

Here,  guided  by  the  results  of  the  previous  trial 
design,  we  change  the  core  properties  to  those  of  1/4  -  5052  -  .002 
in  Figure  5.  The  optimized  design  for  these  properties  is 
shown  on  Figure  6h. 

(i)  Figure  6i 

This  figure  contains  the  data  changes  and  the 
corresponding  results  for  core  1/4  -  5052  -  .0015. 

( j )  Figure  6j 

This  figure  contains  the  data  changes  and  the 
corresponding  results  for  core  1/4  -  5052  -  .0025. 

(k)  Figure  6k 

The  core  cell  size  and  foil  thickness  computed 
(.240961"  and  .002705")  in  Figure  6j  are  very  near  the 
corresponding  values  for  the  1/4  -  5052  -  .0025  core  whose 
material  properties  were  input  for  that  trial  design. 

Therefore,  we  accept  the  1/4  -  5052  -  .0025  core  as  the  final 
core  to  use  in  the  optimum  design. 

Figure  6k  illustrates  the  use  of  the  short  form  of  the 
data  input  in  which  we  specify  only  three  design  parameters,  the 
two  face  thicknesses  and  the  core  thickness.  The  cell  size  and 
the  foil  thickness  are  preset  to  .25"  and  .0025",  respectively. 


All  other  data  are  the  same  as  for  the  Figure  6j  trial  design. 

The  final  optimized  panel  design  parameters  are  shown  on 
Figure  6k.  In  particular,  the  face  thicknesses  should  be  .021904" 
and  the  core  thickness  should  be  1.050198"  to  give  an  optimized 
weight  of  22.980527  lb. 

(1)  Figure  61 

This  figure  illustrates  how  to  exit  from  the  SANOPT 
program,  and  the  logout  procedure. 


SECTION  4 

CONCLUSIONS  AND  RECOMMENDATIONS 


An  approach  to  the  optimum  design  of  flat,  rectangular 
sandwich  panels  subjected  to  edgewise  compression  loading  has 
been  developed.  Using  a  nonlinear  programming  method  called  the 
Sequential  Unconstrained  Minimization  Technique,  the  design 
parameters  characteristic  of  sandwich  construction  are  determined 
so  that  the  final  panel  configuration  is  not  just  an  acceptable 
design,  but  is  the  best  acceptable  design  as  far  as  total  panel 
weight  is  concerned.  The  design  parameters  are  the  individual 
thicknesses  of  the  top  and  bottom  face  sheets  and  of  the  core, 
the  cell  size  and  foil  thickness  of  honeycomb  core,  and  the 
density  of  foam  core.  The  acceptable  ranges  of  the  design 
parameters  are  limited  by  constraints  which  ensure,  for 
example,  that  the  thicknesses  are  within  minimum  and  maximum 
limits,  that  the  faces  do  not  yield,  wrinkle,  or  dimple,  and 
that  the  panel  does  not  buckle.  Possible  panel  loads  are 
edgewise  compression  in  two  directions. 

The  computer  program  which  implements  the  minimum  weight 
sandwich  panel  design  procedure  is  an  easy  to  use  interactive 
code.  Data  are  entered  by  the  user  in  response  to  self- 
explanatory  prompts  from  the  program.  The  user  is  offered 
the  opportunity  to  change  data  before  the  optimization  procedure 
is  initiated.  The  results  of  the  minimum  weight  computations 
are  presented  to  the  user  within  a  few  seconds  in  a  concise 
yet  informative  format.  Typical  output  data  include  the 
optimum  values  of  the  design  variables,  the  minimized  panel 
weight,  and  final  values  of  the  design  variable  constraints. 

The  results  obtained  here  indicate  that  the  minimum 
weight  design  of  sandwich  panels  using  an  analytical  optimization 
procedure  can  be  effective.  The  particular  panel  geometry  and 
loading  considered  are  of  limited  scope  so  that  the  analyses 
required  during  the  optimization  process  can  be  performed 
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using  closed  form  approximate  solutions.  If  the  analytical 
optimization  concept  is  coupled  with  an  efficient  finite 
element  analysis  procedure,  then  more  complex  structural 
configurations  can  be  considered. 

A  useful  extension  to  the  work  presented  here  would  be 
to  integrate  the  sandwich  panel  optimization  procedure  into 

g 

a  multilevel  design  approach.  The  multilevel  design  technique 

q 

combines  the  advantages  of  optimality  criteria  approaches  and 
of  constrained  optimization  procedures  such  as  the  one  presented 
here  for  the  design  of  sandwich  panels.  If  the  multilevel 
design  philosophy  were  adopted,  then  for  example,  complete 
aircraft  wings  could  be  designed  for  minimum  weight.  In  this 
case  sandwich  panels  would  be  only  a  part  of  the  total  structure. 
It  is  recommended  that  an  effort  be  undertaken  to  extend  the 
existing  methodology  for  optimization  of  structural  sandwich 
designs  to  more  general  structures  including  sandwich  having 
laminated  composite  faces,  and  aircraft  wing  structures 
utilizing  sandwich  components. 
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